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ABSTRACT 

The dynamic evolution of the baryonic intergalactic medium (IGM) caused by the underlying 
dark matter gravity is governed by the Navier-Stokes equations in which many cooling and 
heating processes are involved. However, it has long been recognized that the growth mode 
dynamics of cosmic matter clustering can be sketched by a random force driven Burgers' equation 
if cooling and heating are ignored. Just how well the dynamics of the IGM can be described 
as a Burgers fluid has not been fully investigated probably because cooling and heating are 
essential for a detailed understanding of the IGM. Using IGM samples produced by a cosmological 
hydrodynamic simulation in which heating and cooling processes are properly accounted for, 
we show that the IGM velocity field in the nonlinear regime shows the features of a Burgers 
fluid, that is, when the Reynolds number is high, the velocity field consists of an ensemble of 
shocks. Consequently, (1) the IGM velocity v is generally smaller than that of dark matter; 
(2) for the smoothed field, the IGM velocity shows tight correlation with dark matter given by 
v ~ svdm, with s < 1, such that the lower the redshift, the smaller s; (3) the velocity PDFs are 
asymmetric between acceleration and deceleration events; (4) the PDF of velocity difference Av = 
v(x + r) — v(x) satisfies the scaling relation for a Burgers fluid, i.e., P(Av) = (l/r y )F(Av/r v ). 
We find the scaling function and parameters for the IGM which are applicable to the entire scale 
range of the samples (0.26 - 8 h" 1 Mpc). These properties show that the similarity mapping 
between the IGM and dark matter is violated on scales much larger than the Jeans length of the 
IGM. 

Subject headings: cosmology: theory - large-scale structure of universe 



1. Introduction 

The evolution of large scale structure in the uni- 
verse is largely governed by the gravitational clus- 
tering of dark matter. The collapsing of the dark 
matter halos form the potential wells that trap 
baryonic gas, or the intergalactic medium (IGM), 
and form light-emitting objects. Clearly under- 
standing the dynamics of the IGM embedded in 
the underlying dark matter field is essential in 
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large scale structure study. 

The cosmic baryonic gas is a Navier-Stokes 
fluid. Unlike an ordinary fluid, the dynamic be- 
havior of the IGM is dominated by the growth 
modes of dark matter. It has long been recog- 
nized that the growth mode dynamics of cosmic 
baryonic gas can be approximately described by 
the random- force-driven Burgers' equation. Zel- 
dovich's pancake theory and the adhesion model 
first indicated that the non-linear evolution of the 
gravitational clustering of cosmic matter could 
be sketched by a Burgers' equation (Gurbatov, 
Saichev, & Shandarin 1989; Vergassola et al. 
1994). By considering cosmic matter as dissi- 
pative, Berera & Fang (1994) showed that the dy- 
namical equations of the velocity fields of cosmic 
matter are essentially a variant of the random- 
force-driven Burgers equation or the KPZ equa- 
tion (Kardar et al 1986). Later, with a two- 



component (dark matter and IGM) generalization 
of the adhesion model, it was found that the veloc- 
ity potential of the baryonic matter is described by 
the Burgers equation driven by the gravitational 
potential of dark matter (Jones, 1999; Matarrese 
and Mohayaee 2002). 

Observations also reveal that the IGM shows 
features that hint at a field governed by the dy- 
namics of Burgers' equation. For example, the 
transmitted flux of QSO's Lya absorption is found 
to be significantly intermittent and its probabil- 
ity distribution functions (PDF) arc long tailed 
and probably lognormal (Jamkhcdkar et al. 2000, 
Pando et al. 2002, Feng et al. 2003, Jamkhedkar 
et al. 2003). This is consistent with the intermit- 
tent spatial and temporal behavior of a developed 
Burgers fluid (Polyakov 1995, Bouchaud, Mezard 
& Parisi 1995, Gurarie & Migdal 1996, Balkovsky 
et al. 1997, Yakhot 1998, Frisch, Bee & Villone 
1999). Moreover, if the random-force in the Burg- 
ers' equation for the baryonic matter arises from 
the gravity of dark matter, the statistical prop- 
erties of the baryonic matter in the non-linear 
regime will decouple from those of the underlying 
dark matter. This mechanism leads to a breaking 
of similarity between the IGM and dark matter 
on scales larger than the Jeans length (He et al. 
2005). This large scale dis-similarity has also been 
found from observations of the entropy floor (Pon- 
man, et al. 1999), X-ray luminosity-temperature 
relation of galactic clusters; the baryon fraction of 
clusters (Ettori & Fabian 1999, Sand et al. 2003); 
and the multi-phase IGM temperature field (Zhan 
et al. 2003) etc. 

However, there is a difference between the the- 
oretical models and the observations. In all the 
above-mentioned theoretical approaches, heating 
and cooling processes of the IGM are ignored. 
That is, the dynamical equation of the growth 
modes of the IGM without heating and cooling can 
be approximated by a Burgers' equation. How- 
ever, heating and cooling are essential for under- 
standing all details of the IGM. Therefore, an im- 
portant question is how well the dynamics of IGM 
can be described as a Burgers fluid that has un- 
dergone significant heating and cooling evolution. 
To investigate this question, we focus on the veloc- 
ity field or velocity potential because predictions 
based on Burgers dynamics directly concern the 
velocity and velocity potential. We use samples 



of the IGM velocity field produced by cosmologi- 
cal hydrodynamical simulations in which all heat- 
ing and cooling processes are properly accounted 
for and investigate whether the IGM displays the 
common features of a Burgers' fluid. This prob- 
lem is important not only for understanding the 
non-linear behavior of baryonic gas, but also for 
the theory of a Burgers' fluid. 

The paper is organized into the following sec- 
tions: §2 reviews the derivation of the Burgers' 
equations and its basic predictions for the bary- 
onic gas. §3 gives a brief description of the hy- 
drodynamic cosmological simulations. §4 gives the 
statistical analysis of the velocity field of the bary- 
onic gas. Finally, in §5, we will present our conclu- 
sions followed by a discussion of the results. The 
velocity field in the discrete wavelet (DWT) rep- 
resentation is given in the Appendix. 

2. Cosmic baryonic gas velocity field 

2.1. Burgers' equation for cosmic baryonic 
gas 

The derivation of Burgers' equation for cosmic 
baryonic gas has been addressed in Feng et al 
(2004), He et al. (2004) and Pando et al. (2004). 
In summary, the baryonic gas is assumed to be an 
ideal fluid satisfying the hydrodynamic equations 

*! + V(l + *)v = (1) 

^ + (v-V)v = ~Vp-V0 (2) 
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where p, v, E andp are, respectively, the mass den- 
sity, peculiar velocity, energy density and pressure 
of the gas. The term K ra< i in Eq.(3) is given by 
radiative heating and cooling per unit volume of 
the baryonic gas. The gravitational potential (f> in 
eqs.(2) and (3) is given by 

V 2 (j) = AnGa 2 p dm S dm , (4) 

where p~dm{t) and 5dm, are respectively, the mean 
mass density and density contrast of the pertur- 
bations of dark matter. Here we assume that the 
gravitational potential <j) is only given by the dark 
matter mass perturbation. 
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To sketch the evolution of gravitational clus- 
tering we only consider the case where all thermal 
processes are approximated by the polytropic re- 
lations p oc pT, T oc pf-\ or T = T (l + 
where <5(x, i) = [p(x, i) — p(t)]/p(t) is the baryon 
mass density perturbation. The momentum equa- 
tion is then 

— + (v-V)v = — i 77— FT-W 5 

at pm p (1 + o) 

where the parameter p is the mean molecular 
weight of the IGM particles, and m p the proton 
mass. Using the linear approximation for 8 in the 
temperature-dependent term of eq.(5), we have 

^ + ( v .V)v = -Vv-V* (6) 

where the coefficient v is given by 
jk B T Q 
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in which D(t) describes the linear growth behav- 
ior. The term with v in eq.(6) describes a dif- 
fusion characterized by the Jeans length k 2 = 
(a 2 /t 2 )(pm p / 1 k B T ). 

To understand cosmic large scale structure for- 
mation, only the growth modes are of interest. In 
this case, the velocity field of the perturbations is 
irrotational and we can define a velocity potential 

by 



v = V<p. 

a 

Substituting eq.(8) into eq.(6), we have 
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Equation (9) is the stochastic-force-driven Burg- 
ers' equations or the KPZ equations (Kardar, 
Parisi & Zhang 1986; Barabasi & Stanley 1995). 
This is the simplest nonlinear Langevin equation 
for growth modes. The second term on the l.h.s. is 
from convection which is the lowest-order term for 
nonlinear growth of clustering. The third on the 
l.h.s. describes relaxation of the clustering by dif- 
fusion. The term on the r.h.s is the gravitational 
potential of dark matter, which provides the ini- 
tially random perturbations. Clustered structures 
will develop from the initial seeds via the compe- 
tition of the convection and diffusion. The Burg- 
ers' equation has been extensively applied in the 



non-linear dynamics of growing structures from a 
smoothed surface driven by random forces (e.g. 
Vicsek 1992). Eq. (9) shows that structure forma- 
tion of baryonic gas in the universe can be thought 
of as the evolution of surface roughening, i.e., an 
initially smoothed distribution evolves into a wrin- 
kled one. 

2.2. Baryonic gas as a Burgers' fluid 

Much theoretical work on Burgers fluids has 
been done in the hydrodynamic community in the 
past decade (Polyakov 1995; Bouchaud et al. 1995; 
Yakhot 1998; Lassig, 2000; Davoudi et al. 2001). 
Here we summarize the salient points necessary 
for testing whether cosmic baryonic gas behaves 
as Burgers' fluid. 

In the Burgers' equation, there are two scales: 
the dissipation length or the Jeans length 1/fcj 
and the correlation length of the random force, 
r c . The intensity of the random force <p can be 
quantified by the density contrast of dark matter 
5dm ■ The basic feature of the Burgers equation is 
that turbulence will developed in the fluid if the 
following condition holds (e.g. Lassig 2000; Feng 
et al 2003) 



(kjr c ) 2 ' 3 (S 2 d J 1/3 > I- 



(10) 



This condition corresponds to the system being in 
the non-linear regime of the dark matter field. In 
linear regime, Sdm <C 1, and therefore, perturba- 
tions on scale r c > 1/kj will not cause Burgers' 
turbulence in the IGM. In the nonlinear regime, 
Sdm > 1, arid Burgers' turbulence on scales larger 
the Jeans length can develop. Since the dark mat- 
ter field becomes nonlinear first on small scales and 
then on large scales, one can expect that Burgers' 
turbulence in IGM develops from small to large 
scales 

Burgers turbulence is qualitatively different 
from turbulence described by the Navier-Stokes 
equations. The latter generally consists of vortices 
on various scales, while the former is a collection 
of shocks. These features arise because for growth 
modes, the fluid is potential and the velocity field 
is irrotational. If v — > 0, the velocity field in the 
Burger prescription acquires singularities due to 
the discontinuities caused by strong shocks. The 
nonlinear feature of the baryonic gas velocity field 
with fully developed Burgers turbulence can be 
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understood as a field consisting of these shocks 
(Lassig 2000). 

Since dark matter is not affected by Burgers 
turbulence, the IGM velocity field dynamically 
decouples from the dark matter field on scales 
larger than the Jeans length when Burgers tur- 
bulence develops. Jeans diffusion will also lead 
to the decoupling between the velocity fields of 
IGM v and dark matter Wdm- However, the de- 
coupling given by the Jeans diffusion is very dif- 
ferent from that given by the Burgers turbulence 
shocks. For the former, the distribution v and 
Av = v(x + r) — v(x — r) is symmetric with re- 
spect to the transformation Vd m —* — Vd TO , i.e., the 
velocity PDF with Av • Vd m /\vdm\ > (accelera- 
tion in the direction of Vdm/\vdm\) is the same as 
the PDF with Av- Vd m /\vdm\ < (deceleration in 
the direction of v d m/Wdm\)- For Burgers turbu- 
lence, the shocks consist of an acceleration ramp 
followed by a rapid deceleration. Therefore, the 
IGM is not symmetric between the sections accel- 
erating and decelerating. Since the acceleration is 
due to the gravity of dark matter, the acceleration 
is generally in the direction of v dm - Thus, the IGM 
velocity field v is asymmetric with respect to the 
transformation v dm — ► — Vd m . The PDF of Av 
with Av • •Vd m /\'Vdm\ > will not be the same as 
the PDF of Av with Av • v dm /\v dm \ < 0. 

If v = 0, Burgers equation (9) is scale free 
when the random force <fr is scale-free. In this 
case, the velocity field of the Burgers fluid is self- 
similar and the PDF of the velocity difference Av 
is scale-invariant. The velocity difference can be 
effectively measured by the longitudinal velocity 
difference defined as Av r = vn(r\,t) — Vu(r2,t), 
where vn = v ■ (ri — T%)/r and r = ri — r2 1 . Av r 
measures the acceleration and deceleration along 
direction ri — Y2)/r. A detailed analysis shows 
that for a self-similar velocity field, the PDF of 
Av r is (e.g. Davoudi et al 2001) 

^(A^) = >(^), (11) 

where parameter y depends on the statistical prop- 
erties of the random force. 

Summarizing, if the IGM can be described as 
a Burgers' fluid, its velocity field should show the 
following properties: 1.) the peculiar velocity of 
the IGM at a given point will generally be lower 
than that of dark matter at the same point, 2.) 



the PDF of Av r • Vd m /Wdm\ > will not be the 
same as the PDF of Av r • •Vd m /Wdm\ < 0, and 3.) 
The PDF of Av r scales. We will test for these 
properties in the following sections. 

3. Cosmological hydrodynamic simula- 
tions 

Burgers equation (9) is obtained under the as- 
sumptions that 1) decaying modes are ignored; 2) 
the cooling and heating of the IGM are ignored 
and replaced by a polytropic equation of state, 
and 3) the diffusion term is given in the linear 
approximation. Moreover, the driving force in the 
Burgers equation for the baryonic gas comes only 
from the dark matter gravity which itself under- 
goes evolution. It is not at all clear that the fea- 
tures of a Burgers fluid occur in the IGM for which 
cooling and heating processes are essential. To an- 
swer this question, we determine whether samples 
produced by a fully cosmological hydrodynamical 
simulation have the properties listed at the end of 
last section. 

3.1. The WENO algorithm 

In order to capture the Burgers fluid fea- 
tures, simulations need to be effective in capturing 
shocks and discontinuities. Hence we do not use 
a Lagrangian approach such as the smoothed par- 
ticle hydrodynamic (SPH) algorithm. The reason 
is the SPH scheme does not handle shocks or dis- 
continuities well because the nature of SPH is to 
smooth the fields (e.g. Borve, Omang, & Trulsen, 
2001; Omang, Borve, & Trulsen 2003.) Conse- 
quently, Burgers fluid features are overlooked. 

We instead take a Eulerian approach. A well- 
known problem with the Eulerian algorithm are 
the unphysical oscillations near a discontinuity. 
An effective method to reduce the spurious os- 
cillations is given by designed limiters, such as 
the Total- Variation Diminishing (TVD) method 
(Hartcn 1987 et al.) or the piecewise parabolic 
method (PPM) (Collella & Woodward 1984). 
However, the accuracy in the TVD method de- 
generates to first order near smooth extrema 
(Godlewski & Raviart 1996). This problem will 
cause errors in calculating the temperature and 
entropy changes because they are determined by 
the difference of the thermal energy P/(j — 1) on 
two sides of the shock. When the Mach number 
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of gas is high, the thermal energy P/("/ — 1) is 
very small compared to the kinetic energy pv 2 /2. 
To overcome this problem, the Essentially Non- 
Oscillatory (ENO) algorithm and Weighted Es- 
sentially Non-Oscillatory (WENO) were proposed 
(Harten et al. 1987; Shu 1998; Fedkiw, Sapiro 
& Shu 2003; Shu 2003). They can simultane- 
ously provide a high order precision for both the 
smooth part of the solution and the sharp shock 
transitions (Liu, Osher, & Chan 1994; Jiang & Shu 
1996). 

WENO has been successfully applied to hydro- 
dynamic problems containing shocks and complex 
structures. WENO has also been used to study 
astrophysical hydrodynamics, including stellar at- 
mospheres (del Zanna, Velli & Londrillo 1994 ), 
high Reynolds number compressible flows with su- 
pernova (Zhang et al. 2003), and high Mach num- 
ber astrophysical jets (Carrillo et al. 2003). In 
the context of cosmological applications, WENO 
has proved to be especially adept at handling the 
Burgers' equation (Shu 1999). Recently, a hybrid 
hydrodynamic/N-body code based on the WENO 
scheme was developed. It has passed typical reli- 
ability tests including the Sedov blast wave and 
the formation of the Zeldovich pancake (Feng, 
Shu & Zhang 2004). This code has successfully 
produced QSO Lya transmitted flux samples, in- 
cluding the high resolution sample HP1700+6416 
(Feng, Pando & Fang 2003). The statistical fea- 
tures of these samples are in good agreement with 
observed features not only on second order mea- 
sures, like the power spectrum, but also up to or- 
ders as high as eighth for the intermittent behav- 
ior. 

3.2. Samples 

For the purpose of this paper, we run the 
hybrid A^-body/hydrodynamic code to trace the 
cosmic evolution of the coupled system of both 
dark matter and baryonic gas in a flat low den- 
sity CDM model (ACDM), which is specified by 
the cosmological parameters (f2 m , ^a, h, cr 8 , fit) = 
(0.3,0.7,0.7,0.9,0.026). The baryon fraction is 
fixed using the constraint from primordial nucle- 
osynthesis as n b = 0.0125/i~ 2 (Olive et al. 1991). 
The value of fib is lower than the current result 
from WMAP. However, since we we need to refer 
our previous results (Feng et. al. 2003, Pando et. 
al. 2004, He et. al. 2004), which are based on the 



simulation samples produced before WMAP, we 
use the older value. The baryon fraction is not a 
parameter that enters into Burgers' equation, and 
therefore Burgers' fluids features do not seriously 
depend on the baryon fraction once the gravita- 
tion of the baryon gas is negligible in comparison 
to that of dark matter. 

The simulations are performed in a periodic, 
cubic box of size L 3 =25 h~ 1 Mpc with a 192 3 grid 
and an equal number of dark matter particles. The 
size of the grid is 25/192 = 33/2 8 = 0.129 h" 1 
Mpc. The thickness of the Burgers' shock is on 
the order of the dissipation length, i.e., the Jeans 
length, which is in the range ~ 0.1 — 0.3 h _1 Mpc 
for redshifts z < 4 (Bi et al. 2003). Therefore, the 
resolution of our simulation is enough to capture 
shocks. For strong shocks, the thickness of the 
shock is less than the dissipation length. How- 
ever, strong shocks are generally associated with 
high temperature regions (i.e. 3> 10 4 K), and the 
local Jeans length is larger. Therefore, the resolu- 
tion is still less than the Jeans length. In Feng et 
al 2003, we used two samples with different reso- 
lution. The major statistical results were found to 
be independent of the resolutions. 

The simulations start at a redshift z = 49 and 
the results are recorded at redshifts z =6.0, 4.0, 
3.0, 2.0, 1.0, 0.5 and 0.0. The time step is chosen 
by the minimum value among the following three 
time scales. The first is from the Courant condi- 
tion. The second time scale is imposed by cosmic 
expansion which requires that Aa/a < 0.02 within 
a single time step. The last time scale comes from 
the requirement that a particle moves not more 
than a fixed fraction of the cell size. 

Atomic processes, including ionization, radia- 
tive cooling and heating in a plasma of hydrogen 
and helium of primordial composition (X = 0.76, 
Y = 0.24) are modeled in the same way as in 
Bahcall & Cen (1992). Processes such as star for- 
mation, feedback due to SN and AGN activities 
have not yet been taken into account. A uniform 
UV-background of ionizing photons is assumed to 
have a power-law spectrum of the form J{v) = 
J 2 i x 10~ 21 (iV^f/)~"crg s _1 cm _2 sr _1 Hz _1 , with 
parameter J 2 i — 1-0 and a = 1. The photoioniz- 
ing flux is suddenly switched on at z > 10 to heat 
the gas and reionize the universe. 

For our work here, we randomly sampled 192 
data points along 500 one-dimensional directions 
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(500 x 192 = 9.6 x 10 4 total points.) At each point, 
the mass density pdm and peculiar velocity com- 
ponent along the line of sight Vdm of dark matter, 
and mass density pi gm , peculiar velocity compo- 
nent along the line of sight v and temperature T 
of the baryonic gas are recorded. At times during 
our analysis it will be necessarry to smooth the 
distributions on scales of 33/23 h _1 Mpc, where 
.7 = 3,4,5,6,7. 

4. Statistical analysis of velocity fields 

4.1. Relation between the IGM and dark 
matter velocity fields 

In the linear growth regime, the IGM velocity 
follows the dark matter velocity point by point, 
i.e., v(x, t) = Vd m (x, t) (Bi, Borner, & Chu 1993; 
Fang et al. 1993; Bi 1993, Nusser 2000.) Burgers' 
turbulence will lead to a deviation of v(x, t) from 
Vdm (x, t) . Figure 1 gives the relation between the 
velocities of the IGM and dark matter at each grid 
point of the samples at redshifts 6, 4, 2, 1, 0.5 
and 0, where the data is smoothed on scale 33/2 7 
h -1 Mpc. As expected, the data points are scat- 
tered around the diagonal line v — Vd m - At high 
redshift z — 6, the scatter is smaller, reflecting 
the fact that the fields are still in the quasi-linear 
regime. At z < 1, the scatter becomes significant. 
The scatter at low redshift cannot be explained 
by the Jeans diffusion or random noise as Figure 
1 shows a remarkable butterfly-like configuration, 
i.e., v mostly lies < v < Vd, n if «dm > 0, while it 
lies > v > Vdm when Vdm < 0. In other words, 
the peculiar velocity of IGM \v\ is generally less 
than \vdm\- As mentioned in §2.2, shocks natu- 
rally yield the butterfly-like configuration because 
shocks always lead to a deceleration of the gas with 
respect to the underlying dark matter. 

Figure 1 also shows an envelope of the scatter 
at about v ~ w<jm/4. This means, the velocity de- 
crease due to shocks generally is not larger than 
a factor of 4. This factor is consistent with the 
shock theory of a polytropic gas (Landau & Lif- 
shitz 1959). The strongest shocks of a polytropic 
gas lead to an decrease of the velocity by a factor 
of ^ (7 — 1 ) / (7 + 1 ) . For these samples the value of 
7 is found to be 5/3 (He at al. 2004). Therefore, 
the velocity decrease should not be larger than a 
factor (7+l)/(7- 1) =4. 

Figure 2 shows the relation beteen v and vj m 



at z = for the density ranges: (1) low density 
Pdm < 1, (2) moderate density 1 < p dm < 5, (3) 
high density 5 < pdm- We see that even in the 
low and moderate density region, the scatter of v 
is still butterfly-like. This means that the velocity 
fields are dominated by shock structures not only 
around massive halos, but also in moderate den- 
sity areas. Burgers turbulence consists of shocks 
regardless of whether the density is high. 

Figure 3 shows again the relation between v and 
Vdm, but v and Vd m arc now smoothed on scales 
d = 33.3 jV h" 1 Mpc. We see that the butterfly 
feature is weaker with increasingly smooth scales. 
However, this does not mean that the deviation of 
v from Vdm is smaller. On the contrary, the dis- 
tribution of v still substantially deviates from the 
diagonal line for smoothing scales from 2 to 8 h _1 
Mpc. This happens because when the smoothing 
scale is larger than the mean separation of strong 
shocks, the strong shocks are localized within each 
area. Consequently, the mean velocity in all re- 
gions is always dominated by these strong shocks. 

This point can be more clearly seen in Figure 4, 
which plots the relation between v and Vdm at red- 
shifts 6, 4, 2, 1, 0.5, and for the field smoothed 
on scales d = 8 h _1 Mpc. Figure 4 shows that 
at redshift 6, the velocity distribution basically 
follows v = Vdm with little scatter. As redshift 
decreases, the distribution of v is always tightly 
correlated with Vdm, but the correlation between 
v and Vdm is not the diagonal line v = Vd, n , but 
v = svdm with the slope s < 1. The parameter s 
decreases with redshift, it is 0.9 at z — 6, and 0.4 
at z=0. Essentially, s is a measure of the effect of 
strong shocks at various redshift eras. Since Burg- 
ers' shocks are related to the formation of cusps 
in the underlying dark matter density field, the 
shock growth is parallel with the development of 
non-linear clustering. Therefore, the deviation of 
s from 1 is bigger for lower redshifts. 

4.2. Relations of velocity-density and velocity- 
temperature 

We next analyze the relation between mass den- 
sity and velocity. Figure 5 shows v vs. pdm, v vs. 
Pigm, v dm vs. p dm and v dm vs. p lgm at each grid 
point of the samples at redshift z = 0. At a given 
Pigm or pdm, the scatter of Vd m is much larger than 
v. 
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More interesting is the qualitative difference 
between the scatter of v and Vdm- The scatter 
of the IGM velocity v has a clear pd m (or pi gm )- 
dcpendcnt envelope: the smaller the mass den- 
sities, the smaller the scatter. This is because 
the IGM flow is potential. In low density ar- 
eas, the gravitational potential does not vary sig- 
nificantly. Bernoulli's equation predicts that the 
greatest possible value of fluid velocity for a po- 
tential flow is proportional to the speed of sound 
(Landau & Lifshitz 1959). Since the speed of 
sound generally increases with gas density, the 
maximum IGM velocity, i.e., the velocity at the 
v — pigm envelope, increases with pi gm . Although 
p igm and pdm arc not the same point-by-point, 
a large pdm on average corresponds to large p 
(Pando et al 2004). Therefore, the relation v vs. 
Pdm shows a similar envelope. On the other hand, 
the dark matter velocity Vd m is not affected by the 
IGM shocks and is not constrained by the speed of 
sound of the IGM. There is no density-dependent 
envelope for Vdm- 

Figure 6 plots the relation between the IGM 
velocity v and temperature T at each grid point 
of the samples at redshift z = 0. The scatter of 
the IGM velocity v with respect to T also has a 
clear envelope: the smaller the temperature, the 
smaller the scatter. This envelope develops for the 
same reason as above: the IGM velocity cannot be 
much larger than the speed sound. However, figure 
6 shows that the envelope increases with T quickly 
at T < 3, 000 K, but slowly at T > 3, 000 K. This 
is probably because most hydrogen atoms are neu- 
tral at T < 3, 000 K and ionized at T > 3, 000 K. 
On the other hand, the velocity of dark matter Vd m 
is not subject to this constrains. Therefore, there 
is no envelope in the Vd m -T scatter distribution. 

It should be pointed out that the envelopes 
shown in Figures 5 and 6 significantly deviate from 
the usual relation v oc \f(T) for the speed of 
sound. The reason is that the baryon gas exists 

in different regions of T p space. That is, the 

baryon gas is in multiple phases. The multiple 
phases of cosmic baryon gas has been addressed 
in He, Feng & Fang (2004) where it is shown that 
the T — p, and S'(entropy) — T and S — p relations 
are significantly multiply phased either at low or 
high temperature or density. The deviation from 
the usual relation of the velocity and temperature 
shown in Figs 7 and 8 is evidence that the baryon 



gas is substantially multiple phased. 

4.3. Asymmetry of the PDF of velocity 
difference 

As mentioned in §2.2, the features of a Burg- 
ers fluid are effectively described by the statistics 
of the longitudinal velocity difference Aw r (x) = 
U||[x + (r/2)] - W|j[x - (r/2)]. This quantity 
describes the velocity difference at physical po- 
sition x and on distance r — |r|. That is, 
Aiv(x) > corresponds to an acceleration event 
in direction r, while Aw r (x) < is a deceleration 
event. More useful is to define the acceleration 
and deceleration with respect to the velocity di- 
rection of the underlying dark matter velocity 
field. For one-dimensional fields, the events of 
Av r (x)vdm(x) /\vdm(x)\ > correspond to an ac- 
celeration with respect to the dark matter velocity 
field, while Av r (x)vd m (x) /\vd m (x)\ < the decel- 
eration events. 

The velocity difference v(x + r) — v(x — r) con- 
tains position x and scale r. It is convenient 
to measure Av r (x) by the discrete wavelet trans- 
form (DWT) Avj t i (see Appendix). Moreover, the 
DWT modes are orthogonal and effectively avoid 
false correlations and scalings. The smoothing can 
also be done with the DWT, as the scaling func- 
tions on various scales are window functions on 
that scale. For our 1-dimensional samples, the 
DWT variable Avjj (see Appendix) describes the 
velocity difference at the physical position of cell 
I and on distance r = 33.3/2-? fr 1 Mpc. With the 
DWT variables, the acceleration and deceleration 
events with respect to the dark matter field are 
given by Av it i{vj t i/\vj^\)d m - 

We showed in §4.1 that the IGM velocity is gen- 
erally lower than the dark matter velocity because 
of shocks. Shocks cause a decrease in the velocity 
in the direction of the gas flow v. This means the 
PDF of Avjd will be asymmetric between acceler- 
ation and deceleration events. Figure 7 plots the 
PDFs of Avjj for samples at z — and at smooth- 
ing scales j — 3, 4, 5 and 6. It is clear that the 
PDFs of acceleration and deceleration branches 
of the distribution are asymmetric. For j < 5, 
the PDFs of deceleration events generally have a 
longer tail than that of acceleration events, i.e., 
there are more big deceleration events than accel- 
eration events. Using a K—S test, we find that the 
probability the the acceleration and deceleration 
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branches are drawn from the same distribution is 
0.048 for j = 4, and 0.28 for j = 3. Therefore, the 
asymmetry is weaker on larger scale. 

Generally, the asymmetric long tail is signifi- 
cant at \Av jt i\ > 20km/s. For | Avji{vjj/\vjj\) dm \ < 
20 km/s, the asymmetry is weaker, because the 
velocity difference caused by the Jeans diffusion 
is equal to about 20-30 km/s. This is also the 
reason that the asymmetry almost disappeared 
on scales equal to and less that 0.5 h" 1 Mpc, or 
j = 6. The K-S probability is 0.f5 for j = 5 and 
0.79 for j = 6. The thickness of shocks due to 
the Jeans diffusion is on order ^/(v 2 ) 1 / 2 , where 
(v 2 ) 1 / 2 is the rms of velocity field. Therefore, 
the velocity difference Av r (x) with r less than 
vjiv 2 ) 1 ! 2 will not be sensitive to the shock dis- 
continuity. Moreover, the 1-dimcnsional samples 
will also weaken the asymmetry. Since Vd m is a 
projection of the 3-dimcnsional Vdrm the asymme- 
try will be missed for events perpendicular to the 
direction of f-dimension sampling. 

Figure 8 shows the PDFs of Avjj on scale j = 
4, but for redshifts z = 0, f . 0,2.0 and 4.0. This 
figure shows that the PDFs is already asymmetric 
at z = 4. The number of large deceleration events 
is higher for lower redshifts. This yields the slope 
s < f in Figure 4 that is smaller for lower redshifts. 

4.4. Scaling of the Av jt i PDFs 

An important property of the Burgers fluid is 
the self-similarity of the velocity field as given by 
eq.(lf ). A simplified version of eq.(ll) can be ob- 
tained as follows. For the IGM field, the veloc- 
ity gradient is (Av r )/r. If the field of the Burg- 
ers' turbulence shocks is self-similar, the PDF of 
(Au r )/r, F[(Av r )/r], should be scale- invariant. In 
this case, the PDF of Av r , P(Av r ), is given by 
P(Av r )d(Av r ) = F[(Av r )/r]d(Av r /r). Thus, we 
have P(Av r ) = (l/r)F[{Av r /r). This is the PDF 
of eq.(ll) with y = 1. 

Using DWT variables and noting that the phys- 
ical scale r is given by index j as r — 33.33/2- 7 , 
eq.(II) can be rewritten as 



P(Av 3 .i) = 2 yj F(2 yj A Vj< i). 



(12) 



Eq.(I2) requires that the dynamics of the velocity 
difference Avj.i be scale-free. It is not immedi- 
ately clear whether the IGM velocity field is really 
scale-invariant because the IGM dynamical equa- 
tions (1) - (4) contain many physical scales, such 



as the Jeans scale; the scales related to cooling 
and heating, etc., and these are not necessarily 
scale free. 

We tried to fit the distribution P(Avj.i) with 
a scale-invariant function like (12). The first re- 
sult is shown in Figure 9, which gives the P(Avjj) 
against Avj i for the deceleration branch for sam- 
ples at z = 0. We found the best scaling function 
is 



P{Av 3 .i) = 0.012 x2°- 73 exp 



2 0.7j 



Aw 



3,1 



50km/s 



0.75 



(13) 

The PDF eq.(13) is consistent with eq.(12). The 
quantity 0.012 is given by the overall normaliza- 
tion and can be absorbed by a unit transforma- 
tion. The scaling behavior of the PDF depends 
on the two scaling factors 2°^ and (2 - 7 -?) - 75 , 
and therefore, the scaling is given by two param- 
eters: y = 0.70, and the index 0.75. The fac- 
tor 50 km/s makes the exponential dimensionless, 
but is not involved in the scaling and has no j- 
dependence. The scaling function (13) excellently 
fits PDF P(Avj.i) on scales j = 2-7, i.e., from the 
highly nonlinear scale 0.26 h _1 Mpc to the weak 
nonlinear scale 8 h" 1 Mpc. 

Figure 10 is similar to Figure 9, but for the 
acceleration branch. The scaling function is now 



2 0.47j 



A't 



0.75 



P(A Vj<l) = 0.028x2°^ ex P y . 

(14) 

Similar to eq.(13), the scaling behavior of eq.(14) 
depends on two parameters: y = 0.47, and the in- 
dex 0.75. The factor 28 km/s makes the exponen- 
tial dimensionless. Eq.(14) also gives a very good 
fit to the PDFs from j = 2 to j = 6. The PDFs 
at j = 7 shows small deviation from eq.(14). This 
is probably because the scale j = 7, or 0.26 h _1 
Mpc is comparable with the Jeans length. The 
parameter y = 0.47 in eq.(14) is significantly dif- 
ferent than value given in Eq.(13). This difference 
indicates that although both the the acceleration 
and deceleration branches scale, the scaling behav- 
ior is not symmetric between the acceleration and 
deceleration branches. 

Both functions (13) and (14) are consistent with 
eq.(12). Burgers' equation predicts that the accel- 
eration and deceleration branches will have differ- 
ent PDFs. All these results show that the scale in- 
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variant behavior of IGM velocity field is the same 
as a Burgers' fluid. According to the general the- 
ory of Burgers' fluids (Davoudi et al 2001), the 
parameter y should be dependent on the statisti- 
cal property of the random force <fr and the rela- 
tion cannot be found analytically. The simulation 
box of size L 3 =25 h _1 Mpc is not much larger 
than 8 h _1 Mpc, but the scaling shown in Figures 
9 and 10 still gives a good fit on scale of 8 h _1 Mpc. 
The reason for this is the use of the velocity dif- 
ference which suppresses the contribution of long 
wave length perterbations. 

5. Discussions and conclusions 

The initial mass density fluctuations of the 
universe are probably Gaussian with a scale-free 
power spectrum. Therefore, in the linear regime 
the dark matter and IGM mass and velocity fields 
are self-similar. In the non-linear regime, the dy- 
namical behavior of the dark matter and IGM 
are very different even though they are coupled 
via gravity. Since the dynamics of the dark mat- 
ter is governed only by the gravitational interac- 
tions which do not have a preferred scale, the dark 
matter should undergo self-similar evolution even 
in the non-linear regime. The evolution of the 
IGM seems to be much more complicated because 
it sensitively depends on thermal, radiative and 
other non-gravitational processes. Consequently, 
the properties of the Burgers' equation of the dark 
matter field cannot simply be transferred to the 
IGM (He et al. 2004.) 

Nevertheless, we showed that the main dynam- 
ical feature of the IGM velocity field can be ex- 
plained as a Burgers fluid. The IGM samples, pro- 
duced by a full hydrodynamic simulation in which 
thermal processes are involved, show all features 
of a Burgers fluid. This is probably because the 
scales considered are larger than the Jeans length, 
but smaller than the scale at which the nonlinear 
evolution is onset. In this range the IGM is non- 
linear, but not sensitive to the scales of thermal 
processes. Therefore, it can also be approximately 
described as a Burgers' fluid. 

However, as a Burgers' fluid, the IGM velocity 
field is not a simple mirror of the underlying dark 
matter velocity field. For instance the PDFs of 
the velocity difference of dark matter is lognormal 
(Yang et al. 2001) and significantly different from 



eqs.(13) and (14). The PDF of the velocity differ- 
ence of dark matter is symmetric between the ac- 
celeration and deceleration branches with respect 
to the velocity field of dark matter, but asymmet- 
ric for the IGM. All these differences come from 
shocks which do not develop in the dark matter. 
Non-linear evolution will inevitably lead to a dy- 
namical break of the similarity between the dark 
matter and gas. 
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A. The DWT variables of velocity field 

For the details of the mathematical properties of the DWT refer to Mallat (1989a, b), Meyer (1992), 
Daubechies (1992), and for physical applications, refer to Fang & Thews (1998). For us here, the most 
important properties are 1.) orthogonality, 2.) completeness, and 3.) locality in both scale (r) and physical 
position (x). Wavelets with compactly supported basis are an excellent means to analyze the velocity 
fields. Among the compactly supported orthogonal wavelets, the Daubechies family of wavelets are easy to 
implement. 

To simplify the notation, we consider a 1-D velocity field v(x) on spatial range L. It is straightforward to 
generalize to 3-D fields. In DWT analysis, the space L is chopped into V segments labeled by I = 0, 1, ...2 J ' — 1. 
Each of the segments has size LjV . The index j is a positive integer which represents scale LjV . The index 
I gives position and corresponds to spatial range ILjV < x < (I + \)Lj2K 

DWT analysis uses two functions, the scaling functions 4>j,i(x) — (2 J /L) 1 / 2 ^^- 5 /L — I), and wavelets 
i>j,l(x) — (2- 7 '/L) 1 / 2 ^>(2- 7 '/L — 0- The scaling functions play the role of window function. They are used to 
calculate the mean field in the segment I. The wavelets i>j,i(x) capture the difference between the mean 
fields at space ranges ILjV < x < (I + l/2)L/%> and (I + l/2)L/2^ < x < (I + \)Lj2K 

The scaling functions and wavelets ipjj(x) satisfy the orthogonal relations 

J <t>j,l ( X )<t>j,l' ( X ) dx = S W i ( Al ) 

J i>iA x )^,v( x )d x = 5 i,i' 5 i,Vi ( A2 ) 
J <t> j ,i{x)i, jlt v{x)dx = 0, if j' > j. (A3) 
With these properties, a 1-D random field v(x) can be decomposed into 

oo 2 j '-l 

v(x) = v j (x) + ej',i*l>j',i(*)> ( A4 ) 

3'=3 1=0 

where 

2 3 -l 

The scaling function coefficient (SFC) ejj and the wavelet function coefficient (WFC), Zjj are given by 

tj,i = J v(x)(t> jt i(x)dx, (A6) 



and 



j = J v{x)i) 0il {x)dx, (A7) 



respectively. The SFC ejj measure the mean of v(x) in the segment I, while the WFC Zjj measures the 
fluctuations (or difference) of field v(x) at I on scale j. 

The first term on the r.h.s. of eq.(A4), i^'(x), is the field v(x) smoothed on the scale j, while the second 
term contains all information on scales > j. Because of the orthogonality, the decomposition between the 
scales of < j (first term) and > j (second term) in eq.(A4) is unambiguous. 

Since scale r, and position x, correspond to, respectively, j, I, and <t>j,i(x') is window filter on scale r at 
position x, the DWT variables of velocity field are given by 

$ v{x)4> 3 .i(x)dx 

V " = f M x)dx ■ (A8) 



10 



Vjj is the mean velocity in the spatial range IL/2 3 < x < (I + l)L/2 3 . Similarly, the DWT counterpart of 
the velocity difference is 

/ v(x)tpj t i(x)dx 
J (j>j,i{x)dx 



Avj t i = J TJ '/Z ■ (A9) 



Avjj is the difference between the mean velocities of spatial ranges /i/2- 7 < x < (I + l/2)L/2 : > and (I + 
l/2)L/2i < x < (I + \)Lj2K 

Vj i and Avj are the variables of the velocity field v(x) in the DWT representation. These variables give 
a complete description of the field v(x) without a loss of information. The orthogonality of scaling functions 
and wavelets insure that the decomposition does not cause false correlation among these variables. 
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Fig. 1. — v vs. Vdm at each grid point of samples velocity difference Av r = v(x + r) — v(x) are taken 

at redshifts 6, 4, 2, 1, 0.5 and 0. The two solid to be j = 2,3,4,5 6 and 7. 

lines are, respectively, v = Vdm, and v = Vdm/^- 

Fig. 2. — v and Vd m at z = for different density 
ranges: (1) (1) p dm < 1 (top), (2) 1 < p dm < 5 
(middle), and (3) pdm > 5 (bottom). 

Fig. 3. — v and Vdm for sample z = 0. The velocity 
field is smoothed on scales d — 33.3/2- 7 h _1 Mpc, 
and j = 2,4,6,8. 

Fig. 4. — v and v dm for samples at z = 6, 4, 2, 
1, 0.5, 0, and smoothed on scale 8 h _1 Mpc. The 
solid line is v — svd m - The v-Vd m relation is fitted 
with line v — svd m (dashed line). The parameter 
s is shown in each panel. 

Fig. 5. — The relation of v vs. pdm (top left); 
v vs. pi gm (top right); v dm vs. p dm (bottom 
left), and Vd m vs. pi gm (bottom right) at each 
grid point of the samples at redshift z = 0. 
The envelop of the v-pd m distribution is given by 
v = [log(p dm /0.001)] h , where b = 2.9 (dotted); 3.0 
(dashed) and 3.1 (solid). 

Fig. 6. — The relation of v vs. temperature T 
(top) and Vdm vs. temperature T (bottom) at 
each grid point of the samples at redshift z = 0. 
The envelop of v-T distribution is given by v = 
[g\og(T/l00)] b , where b = 2.5; g = 2.6 (solid), 3.0 
(dashed) and 3.5 (solid). 

Fig. 7. — The PDF of velocity difference Avjj for 
acceleration events (solid) and deceleration events 
(dotted). The samples are at z = with smooth- 
ing scale d — 33. 3/2^ h _1 and j — 3, 4, 5 and 6. 

Fig. 8. — The PDF of velocity difference Avjj for 
acceleration events (solid) and deceleration events 
(dotted). The samples are at redshifts z = 4, 2, 1 
and 0, with smooth scale is j = 4. 

Fig. 9. — The PDF of velocity difference |Auj,j| of 
deceleration cases, i.e. Avjd{vj t i/\vj t i\)d m < The 
samples is for z = 0. The scales r in the velocity 
difference Av r = v(x + r) — v(x) are taken to be 
j = 2,3,4,5 6 and 7. 

Fig. 10. — The PDF of velocity difference \Avjj\ 
of acceleration cases, i.e. Avj t i(vj t i/\vjj\)d m > 
The samples is for z = 0. The scales r in the 



14 



This figure "Fig01.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



This figure "Fig02.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



This figure "Fig03.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



This figure "Fig04.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



This figure "Fig05.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



This figure "Fig06.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



This figure "Fig07.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



This figure "Fig08.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



This figure "Fig09.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



This figure "Figl0.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0502175v2 



